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Abstract. We present results from the XMM-Newton observation of the binary cluster A1750 a,t z — 0.086. We 
have performed a detailed study of the surface brightness, temperature and entropy distribution and confirm that 
the two main clusters of the system (A1750 N and A1750 C) have just started to interact. From the temperature 
distribution, we calculate that they are likely to merge sometime in the next 1 Gyr. The more massive cluster, 
A1750 C, displays a more complicated temperature structure than expected. We detect a hot region associated 
with a density jump ~ 450 kpc east of the cluster centre, which appears to be a shock wave. This shock is not 
connected with the binary merger, but it is intrinsic to A1750 C itself. From simple physical arguments and 
comparison with numerical simulations, we argue that this shock is related to a merging event that A1750 C has 
suffered in the past 1-2 Gyr. The larger scale structure around A1750 suggests that the system belongs to a rich 
supercluster, which would presumably increase the likelihood of merger events. These new XMM-Newton data 
thus show that A1750 is a complex system, where two clusters are starting to interact before having re-established 
equilibrium after a previous merger. This merger within a merger indicates that the present day morphology of 
clusters may depend not only on on-going interactions or the last major merging event, but also on the more 
ancient merger history, especially in dense environments. 

Key words, galaxies: clusters: individual: A1750 - galaxies:clusters: intergalactic medium, mergers - X- 
ray:galaxies:clusters - X-rays:general - 



1. Introduction 

Clusters of galaxies are the largest gravitationally bound 
structures in the Universe, and because of their observable 
properties and sheer number, they represent invaluable 
cosmological probes. 

The hierarchical scenario of structure formation in 
the Universe predicts the growth of clusters of galaxies 
through accretion of smaller units by gravitational infall 
and mergers. During a merger event more than 10^^ ergs 
are dissipated in the intra-cluster medium (ICM) by shock 
heating, leading to strong variations of physical charac- 
teristics such as the temperature, density and entropy. 
Numerical N-body simulations have shown that mergers 
produce substructure which is detected in both the density 
and temperature distributions (Schindler & Miiller [T993I 
Roettiger et al ll997l Ricker ni)98|l . The temperature sub- 
structures survive for ~ 4-6 times longer than density sub- 
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structures, thus the temperature distribution is a strong 
indicator of the cluster history and dynamical state. X- 
ray spectro-imaging observations are well suited to inves- 
tigate in a deeper detail the observable characteristics of 
the merger event. 

Density substructure in X-ray observations of some 
low-redshift clusters (Forman & Jones 119821 Mohr, 
Fabricant & Geller [Tl)93|l . together with the detection of 
other substructure in optical observations (Escalera et al. 
I1994f) , suggests that these clusters are dynamically young 
(Buote & TsaiUninj 

The first quantitative X-ray temperature maps were 
obtained with ASCA (e.g., Markevitch et al ll998|l but, as 
with the later BeppoSax observations (e.g., De Grandi & 
Molendi I1999|l . resolution was low due to the energy de- 
pendent Point Spread Function (PSF) of these telescopes. 
Thanks to the high effective area and spatial resolution 
now available from XMM-Newton and Chandra, precise 
spatially resolved temperature maps are now possible. 
Results confirm the increase in temperature in the merging 
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regions (e.g., Arnaud et al. 121)011 Markevitch & Vihklinin 
IMin Markevitch et al. Neumann et al. 2003 ) but 

the new observations have also detected new phenomena, 
such as the cold fronts (Markevitch et al. l2000|l . 

Substructures in clusters contain a fossil record of the 
merger history. Statistical studies of cluster morphology 
can thus provide an important test of cosmological mod- 
els of structure formation. In practise, attempts to com- 
pare the prediction of numerical simulations of structure 
formation to quantitative statistical studies of nearby clus- 
ter morphology have been hampered by our current poor 
knowledge of the effect of cluster growth on morphology. 
Despite the large amount of data available, the dynamical 
evolution of the ICM and the relation between the galaxies 
and gas during a merger event is still poorly understood 
(see Buote 12002 for a review). A better understanding of 
the physics of merger events, in particular the relaxation 
time of substructures fNakamura ll99^ is required. 

Hence a systematic study of merging clusters of galax- 
ies would help a great deal in deepening our understand- 
ing of the process of cluster formation and evolution. Our 
XMM-Newton Guaranteed Time (GT) time program was 
established for this purpose. We selected a small sample of 
clusters displaying the signature of substructure (on the 
basis of ROSAT imaging analysis), and which could be 
considered as being in different epochs of a merger event. 

In this paper we present the first object in our sample: 
the binary cluster A1750. In optical, A1750 shows a multi- 
peaked galaxy distribution (Beers et al. 119911 Ramirez & 
Quintana HMH Donnelly et al. lM)T|l : the two major peaks 
have a radial velocity difference of ~ 1300 km s~^. A1750 
was observed in X-ray with Einstein (Forman et al. 11981(1 . 
and subsequently identified as a canonical binary cluster, 
though three peaks are clearly visible in the Einstein X-ray 
image. The cluster was later observed with ROSAT and 
ASCA (Novicki et al. ITMHl Donnelly et al. EOOB; three 
peaks are also visible in the 2° ROSAT image. The latter 
authors combine ROSAT and ASGA observations, as well 
as the galaxy distribution, giving the first kinematic and 
dynamical description of the cluster. They detected an 
enhancement of order 30% in the temperature between the 
two main peaks in the X-ray emission, and suggested that 
a shock region is developing in the gas by compression. 
However, Donnelly et al. H2001|l concluded that the two 
clusters could effectively be considered as isolated objects, 
as the hot gas in the interacting region did not have a 
strong impact in the global temperature estimate, and its 
contribution to the total density was weak. Donnelly et 
al.'s mass analysis indicated a mass ratio of order 1.3, but 
their dynamical analysis was not conclusive as to whether 
or not the system is bound. 

In this paper we use the high sensitivity and spatial 
resolution of XMM-Newton to gain new insights into this 
system. Taking advantage of the large field of view of 
XMM-Newton^ we observed the two main subclusters in 
a single pointing. The exceptional sensitivity of XMM- 
Newton allows us to produce temperature and entropy 
maps with an accuracy which considerably surpasses pre- 



vious attempts. We detect new features in all maps and 
relate these to the dynamical state of the system, arguing 
that, while the clusters are just beginning to interact, at 
least one of them (A1750 C), and possibly both, may not 
yet have relaxed from a previous merger event. 

Throughout the paper we assume Hq = 50 km 
Mpc~^, q=0.5. In this cosmology 1 arcmin corresponds to 
130 kpc at the mean redshift of the cluster (z — 0.086). 

2. Observations and data preparation 

2.1. Observations 

A1750 was observed for 34 ks in July 2001 (Revolution 
300) by XMM-Newton. In this paper only data from the 
European Photon Imaging Gamera (EPIG; Striider et al. 
120011 Turner et al. I2()()l|l are considered. Galibrated event 
list files were provided by the XMM-Newton SOC. The 
observations were obtained with the MEDIUM filter; the 
full frame mode was used for MOS and the extended full 
frame for pn. Throughout this analysis single pixel events 
for the pn data (PATTERN 0) were selected, while for the 
MOS data sets the PATTERNS 0-12 were used. 

The background estimates were obtained using a 
blank-sky observation consisting of several high-latitude 
pointings with sources removed (Lumb et al. The 
blank-sky background events were selected using the same 
selection criteria (such as PATTERN, FLAG, etc.) used 
for the observation events. Furthermore, the blank-sky 
background file was recast in order to have the same sky 
coordinates as A1750, ensuring that the source and back- 
ground products come from the same region of the de- 
tector, reducing errors induced by any detector position 
dependence. 

The source and background events were corrected 
for vignetting using the weighted vignetting method de- 
scribed in Arnaud et al. (^2001J. This allows us to use the 
on-axis response matrices and effective areas. 

2.2. Background estimate 

The estimate of the background level is a crucial point 
since we are interested in extended and low surface- 
brightness sources, were the background effects are im- 
portant especially at large distance from the centre. 

The XMM-Newton background consists of several 
components, which may be variable in time and in space 
distribution. The soft proton background is a time and 
flux variable component. Within some observations it is 
possible to find several periods of time where the mean 
flux level varies by a factor of order 10 or more. To clean 
the data for this emission we reject all frames outside 2a 
of the mean value in the 10-12 keV (12-14 keV for pn) 
light curves for each camera, using the method of Pratt & 
Arnaud (.2002;) to determine the mean and a values. For 
this analysis, light curves were grouped in bins of 100 s. 
The observation was relatively clean: using these criteria 
~ 10% of the total exposure time was rejected. The final 
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(alb) 

Fig. 1. (a) EPIC counts image in the energy band 0.3-7.0 keV; the image is not corrected for vignetting (b) Iso- 
intensity contours of the M0S1-I-M0S2 image in the energy band [0.3-1.4] keV overlaid on the optical DSS image. 
Contours are logarithmic, and were obtained by generating contours after replacing point sources with Poisson noise 
from surrounding annuli. 



exposure times are 30.4 ks, 31.4 ks, 23.4 ks for MOSl, 
M0S2 and pn cameras respectively. The blank-sky back- 
ground event files were then cleaned using the same crite- 
ria. 



A second component of the XMM-Newton background 
is represented by the particle induced background, which 
dominates at high energy (> 5 keV) and induces fluo- 
rescence lines (Al, Si, Cu, Au) from the shielding of the 
camera and the detector itself. This background compo- 
nent is effectively described by the blank-sky background, 
under the hypothesis that the variation of this particle 
background is small. It has been demonstrated that the 
particle background is variable at the ~ 10%-level. This 
variation is sufficiently small that it can be taken into ac- 
count by a normalisation in the [10-12] keV (MOS) and 
[12-14] keV band (pn) between the source and the blank- 
sky background observations. The normalisation used in 
this work are 1.09, 1.11, 1.11 for MOSl, M0S2 and pn 
respectively. 



Finally, we have to consider the astrophysical X-ray 
background, which is the combination of a soft (E< 1.5 
keV) component, mainly due to the local bubble, and a 
hard component due to unresolved cosmological sources 
(mainly AGN). This emission is properly an X-ray com- 
ponent (i.e., not particle- induced, or due to soft proton 
flares). The hard component of this astrophysical X-ray 
background is well taken into account by the subtraction 
of the blank-sky background. At lower energies, however, 
the astrophysical X-ray background is variable across the 
sky (e.g., Snowden et al. 11997 1. To account for this, one 
possibility is to use a local background as described in 
Pratt, Arnaud, Aghanim ( 2001) and Arnaud et al. (|2002() . 
However, emission from A1750 fills the whole field of view 
(FoV) and there is no large region which can be consid- 
ered free of cluster emission with enough confidence to be 
used as a local background. In this analysis, the blank-sky 
background represents the best estimate we can give of 
the background associated with the cluster. However, the 
examination of radial profiles in several different bands 
suggests that at large radii the emission remaining after 
the subtraction of the blank-sky background is essentially 
of order zero, giving us confidence in our use of the blank- 
sky background only. 



3. Morphology 

3.1. Image 



Figure 1(a) shows the combined MOSl, M0S2 and pn 
image of A1750 in the energy band 0.3-7.0 keV. The im- 
age is in counts, has not been corrected for exposure, and 
the regions outside the field of view (FoV) have not been 
masked. A total of 356000 photons have been collected 
by EPIC in this energy band. The double distribution of 
the emission is clear, with one emission peak at the cen- 
tre of the FoV and the other to the north-east. There are 
also many serendipitous point sources and some extended 
sources: their positions and fluxes were obtained using 
the EMSRLI file in the pipeproducts, and all those with 
flux greater than 10~^^ ergs arcmin"^ were masked 
throughout this analysis. 

Contours of the low-energy (0.3-1.4 keV), adaptively- 
smoothed EPIC/MOS image are shown overlaid on the 
DSS image in Fig. 1(b) The X-ray emission peaks are 
at [a = 13''30™49".881, S = -01°51'46'.'70] and [a = 
13''31"10".941, S = -01°43'41'.'65] for the cluster at the 
centre of the FoV (hereafter A1750 C) and the north- 
eastern cluster (hereafter A1750 N) respectively. 

The isophotes at the centre of A1750 C are elongated 
and their ellipticity decreases with the distance from the 
centre. The distance between the peak of X-ray emission 
of A1750 C and the cD galaxy visible in the DSS image 
is ~ 15" (32.5 kpc). This shift is larger than the possible 
attitude error of XMM-Newton. The cD galaxy is offset 
from the X-ray peak towards the east, in which direction 
there is a clear compression of the isophotes. 

The centroid of the X-ray emission of A1750 N is lo- 
cated exactly over the superimposed optical emission of 
two central galaxies. The X-ray isophotes in the central re- 
gions are quite circular but they show a slight compression 
to the east-south-east. The abrupt change in the shape and 
orientation of the isophotes at large scale is however re- 
markable, where they become very elliptical and oriented 
in the direction of the line joining the centres of the two 
clusters. 



3.2. 2D P model fitting 

To take into account the binary nature of the cluster and 
to give a quantitative measure of the morphological char- 
acteristics we observe, and in an attempt to detect density 
substructures, we performed a bi-dimensional analysis. 
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Table 1. Best fit results of the 2D (3 model. 





A1750 N 


A1750 C 


Tci (kpc) 


250 


230 


rc2 (kpc) 


180 


200 


/3 


0.56 


0.68 


ao (J2000) 


13''31'"10''5 


13''30"'49'0 


So (J2000) 


-01°43f20'.'61 


-01°51'54'.'14 


PA 


-12° 12' 53" 


24°38'14" 



Fig. 2. EPIC/MOS adaptively smoothed image with, 
overlaid, contours of the residuals after the subtraction 
of a double f3 model. The residuals are traced at 2,3,5 and 
7 a above the background. Further details in the text. 



We adopted the hypothesis that a /? model (Cavaliere 
& Fusco-Femiano I1976|l is an appropriate description of 
a relaxed isothermal cluster. We fit the surface brightness 
distribution of A1750 at low energy (where the density dis- 
tribution is least temperature dependent) with a /3 model 
and quantify the deviation from this model and the pres- 
ence of substructures. In fact, the interest of this exercise 
is to compare A1750 with the surface brightness distribu- 
tion of two relaxed clusters. 

For the fit we followed the prescription of Neumann 
& Bohringer (1997) and apply this analysis to the MOS 
camera only because the large gaps in the pn camera are 
relatively difficult to take into account. Images of pixel size 
of 4'.'1 from the weighted MOS event files and in the energy 
band 0.3-1.4 keV were created, and then summed. Since 
the weight correction is applied directly to the events, the 
statistical Poisson distribution of photons does not apply 
for these images. To correctly take into account the errors, 
the variance is given by — 'Siwf where w^ is the weight 
associated to the iih photon (see Arnaud et al. 17001 ). 
The error images were generated for each MOS image and 
the quadratic sum of the latter gives the summed error 
image. The combined MOS image was smoothed with a 
Gauss filter of = 15", chosen to be of the same order of 
the MOS PSF. This smoothing allows us to be confident 
that in each spatial interval of integration (meta-pixel) 
we can use the Gaussian statistic and find the best 
when fitting. The error image was treated according to the 
error propagation function for Gaussian filtered images 
described in Neumann & Bohringer (1997). 

Because of the double nature of A1750, the 2D /3 model 
for a single cluster (see Neumann & B6hringer il997.. Pratt 
& Arnaud l2002|l has been modified in order to take into 
account simultaneously the two clusters. The background 
was approximated as a unique constant value. The region 
between the two clusters was masked before fitting so as to 
reduce the interference from a possible interaction region. 
Point sources were also masked. Results of the best fit dou- 
ble 2D /3-model are shown in Tabled The best fit model 
was then subtracted from the image in order to quantify 
the significance of any possible excess flux. The signifi- 



cance was calculated following the Appendix of Neumann 
& Bohringer (1997). The positive residuals, traced at 2, 3, 
5 and 7 a above the background are shown in Fig. [3 

The image clearly shows significant residuals in the 
centre of each cluster. The core of A1750 N is off-centre 
with respect to the structure at larger scale (which is well 
described by the /3-model). The same effect is marginally 
visible at the centre of A1750 C. Other significant residu- 
als, other than the two cores, are detected to the north of 
A1750 N and to the west of A1750 C. 

This analysis does not show any significant substruc- 
ture in the region between the two clusters (hereafter the 
middle region), with the exception of two marginally de- 
tected (~ 2a) point sources and a further source which is 
detected at greater than Scr. From comparison with the 
raw image, it is clear that this source is extended. 

We have considered the nature of the extended source 
in the middle region. It is clearly detected in emission in 



the pn image (see Fig. 1(a) I, but unfortunately most of the 
emission is located in the gaps of the MOS images. We ex- 
tracted the spectrum of this source in a circle of 45". As 
a background we used an annular region centred on the 
source, with an inner radius 1.'3 and outer radius 2.'2. The 
spectrum does not contain emission lines strong enough to 
estimate the source redshift. We fitted the spectrum with 
an absorbed MEKAL model with redshift and column 
density parameters fixed to those of A1750 (z = 0.086 and 
UH = 2.37 X 10^0 cm-2). The best fit gives kT = i.GtH 
keV and a flux of 1.3 x 10~^^ erg cm^^ s^^ in the energy 
range [0.3-10.0] keV. A superposition of the residual con- 
tours on the DSS image does not show any obvious optical 
counterpart. 



4. Thermal structure 

The lack of significant substructure in the middle region 
seems to suggest that there is little or no interaction be- 
tween the two components of the cluster (however, see 
Sect. I7.l|l . In other words, the merger event is sufficiently 
early that the density enhancement expected in the inter- 
action region (e.g. Roettiger et al. 1997; Ritchie & Thomas 
2002) has not yet been produced (or is only just begin- 
ning). Knowledge of temperature structure is thus a fur- 
ther important step toward the comprehension of the dy- 
namical state of the system. 

The spectral analysis was performed by extract- 
ing spectra from the weighted source and background 
(blank-sky) event lists. The out of time events affect- 
ing the pn observation were not taken into account, 
but in our experience they are not expected to have 
a detectable effect on the results. The response files 
used in this analysis are: ml_mediumv9q20t5r6_alL15.rsp 
(MOSl), m2.mediumv9q20t5r6.all.l5.rsp (M0S2) and 
epn_ef20_sY9 .medium. rsp (pn). 
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Fig. 3. EPIC/MOS adaptively smoothed image in the en- 
ergy range 0.3-0.6 keV with superimposed contours of the 
3 regions of interest for the global spectrum analysis. The 
image is logarithmically scaled. 



4.1. Global spectrum 

We obtained a global spectrum for each of the regions dis- 
played in Fig. 131 The regions were chosen to give a first 
measure of the spectral characteristics of the two clusters 
and the middle region, and also serve as a useful com- 
parison with previous results. The particular shape of the 
middle extraction region is determined by the exclusion of 
point sources and cluster emission. The spectrum of the 
background was extracted in the same region and, after 
normalisation, subtracted from the source spectrum. 

The spectra of the 3 cameras, excluding the channels 
below 0.3 keV, were simultaneously fitted with an ab- 
sorbed MEKAL model with the temperature, metallicity 
and normalisation (emission measure) as free parameters. 
The column density was fixed to the best fit value ob- 
tained by fitting the MOS cameras only. This value is 
in agreement with the galactic value for A1750 C. (The 
spectrum of A1750 N and the middle region are better 
fitted with a lower than Galactic absorption). The best 
fit results are listed in Table |21 The temperature in the 
middle region is significantly higher (by ~ 30%) than that 
of A1750 C, in agreement with the results of Donnelly 
et al. (^001). In order to be confident in this detection 
of a significantly higher temperature, we also investigated 
the dependence of the temperature of middle region spec- 
trum on the Galactic absorption. If we fix the Nh to the 
Galactic value, the temperature drops to 4.35±0.45 keV 
(for a reduced x^=425. 86/427 d.o.f.), while if we fix the 
Nh to the value obtained by fitting cluster N, the best 
fit temperature is 4.9tEJj keV (x^=418.98/427 d.o.f.). In 
both cases the temperature is higher than the temperature 
of A1750 C at the 90% confidence level, and the goodness 
of the fit is clearly in favour of a lower absorption. In Table 
12 we list the temperatures obtained by fitting simultane- 
ously the three cameras with the column density fixed to 
the best fit MOS value in each region. The spectral fits 
are shown in Fig. 0] 

For A1750 C and the middle region, our results are 
in excellent agreement with the ASCA analysis (Donnelly 
et al. 2001), but we find a lower temperature than for 
A1750 N, if the mean temperature between zones 1 and 
2 in Donnelly et al.'s paper is taken. However, if we com- 
pare our result with the ASCA temperature in region 1 
of Donnelly et al.'s Figure 2, the two results are in better 
agreement. 

We further fit the global spectrum of A1750 N with 
a two temperature model. In this case we find kTi = 
3.17 ±0.1 keV, kT2 0.65 ± 0.17, with abundances of 
0.3±0.7 solar and a x^=:754. 05/660 d.o.f. According to the 
F-test, it is 99% probable that the 2 temperature model 
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h-H — j \ \ \ \ \ '■ 



channel energy (keV) 

Fig. 4. Global spectra of the three regions in Fig.|niAl750 
C (Top); A1750 N (centre); middle region (bottom). MOS 
spectra are in black, the pn in grey. The three camera 
spectra were fitted simultaneously. 

is a better representation of the data, suggesting that the 
gas is multi-temperature (see also Sect. El . 

4.2. Temperature map 

4.3. Adaptive binning 

The high sensitivity of XMM-Newton allows us to build 
a temperature map. To optimise the spatial resolution, 
without any spatial a priori, we have defined a specific 
binning which follows the statistics, and leads to fairly 
homogeneous errors in the temperature determination for 
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Table 2. Results of the global spectrum analysis in the region shown in figure |3| The column density is obtained by 
fitting the MOS spectra only and then fixed to this value when combining with the pn spectrum. The value of the 
column density listed here is the best fit MOS value (see text for more details). Errors are at 90% for one significant 
parameter. 



region 


Nh 

( 10^" at/cni-2) 


kT 
(keV) 


Z 
solar 


xVd.o.f 


Lx [0.2-10] keV 
(erg s-i) 


Cluster C 


2.18±0.44 


3.87±0.10 


0.32±0.04 


1217.79/1139 


2.2 10"* 


Cluster N 


1.00±0.8 


2.84±0.12 


0.22±0.05 


797.33/663 


1.3 10"'' 


Middle 


0.5(< 1.4) 


r 1 2+0.77 


n 97+0.28 
•J-^' -0.19 


418.32/426 


1.7 10"^ 



each meta-pixel. This approach is similar to the one pre- 
sented by Sanders & Fabian 1)2001(1 . We produce a full 
energy band image with sources excised. The algorithm 
starts from the largest spatial scale and continues down to 
the smallest possible, dividing by 4 each cell of the image 
containing more than a given number of photons (here we 
choose 2000 photons for the 3 cameras in the full energy 
band prior to background subtraction). We end up with 
a list of cells that have more or less the same number of 
photons. We apply this spatial binning to the observation 
and to the background to extract all spectra. Then using 
XSPEC, we fit the MEKAL model (with the Nh fixed to 
the galactic value and abundances set to 0.3 solar) to each 
spectrum. The fit temperature of each cell is then used to 



tested on simulated observations; here we present prelim- 
inary results from the first application to a real observa- 
tion. For a detailed description of the method sec Bourdin 
et al. 



build the temperature map shown in Fig. 5(a) The typ- 
ical error (which depends slightly on the temperature) is 
around ±0.6 keV. 

As already suggested by the broad global spectrum 
analysis above, the region between the two clusters dis- 
plays a higher temperature, implying a dynamical inter- 
action. The map is in excellent agreement with that of 
Donnelly et al. fSTOUj but because of the improved sen- 
sitivity, our map indicates that the temperature structure 
is more complicated. 

A1750 N displays a relatively uniform temperature 
with a suggestion of a cooler core region (although the 
large bins are not optimised). To the west side of A1750 
N we observe a higher temperature which may be con- 
nected with the interaction between the two clusters. A 
hot region to the south of the centre of A1750 N should 
be also noted, which interestingly lies ahead of the com- 
pression of the isophotes near the core. 

A new feature, observed for the first time in these data, 
is the hot region '^3-4 arcmin to the east of the centre of 
A1750 C. This hot region corresponds exactly to a com- 
pression of the isophotes in the X-ray surface brightness, 
and does not appear to be related to the larger hot re- 
gion due to the interaction between the clusters. In other 
words, it appears to be intrinsic to A1750 C itself. 

4.4. Wavelet temperature analysis 

We also computed a temperature map by applying the 
new multi-scale spectro-imagery algorithm described in 
Bourdin et al. lj5003,l . In that paper, the algorithm was 



The A1750 temperature map was computed using only 
the event lists from the MOS cameras, as we do not yet 
have a stable model for the pn background. Due to the 
steep decrease of the cluster emissivity with radius, the 
analysis has been restricted to regions where significant 
signal is detected at more than one scale. The background 
normalisation was fixed to a value estimated for each cam- 
era in an external region of the field of view, where the 
background emissivity dominates. The temperature map 
is a 128 X 128 pixels image, where the multi-scale analysis 
has been performed on 5 scales in order to prevent large 
scale edge effects, so that details on spatial scales from 
around 15" to 4' are expected. Spectral fits with error-bars 
higher than 1.5 keV are not taken into account due to ev- 
idence for a bias of the temperature estimator in regions 
with high background level . The issue of background mod- 
elling for spectral fits in each element is the major source 
of uncertainty, especially at large distance from the cen- 
tre where the background itself begins to dominate. To be 
conservative, we have applied a posteriori a mask to the 
wavelet maps defined such that at least 10 photons are 
detected in the emissivity wavelet map. This corresponds 
to a 3ct signal above the background. 

Since we have not yet perfected the removal of point 
sources in the algorithm, the following wavelet- derived 
maps are presented without point sources removed. The 
side by side comparison of the adaptive binning and 
wavelet temperature maps (Figs 5a and 5b) allows the 
interested reader to judge the effect of the point sources 
on the temperature structure. 

Note that, in the wavelet map, all the point sources 
appear in the highest frequency plane. This means that 
their contribution to the local temperature is limited to 
that spatial frequency (in other words they always appear 
in the temperature map at the highest frequency, i.e., as 
point-like sources). 

In general, the wavelet-derived temperature (Fig. 5b; 
errors Fig. 5d) map is in good agreement with the 
adaptively-binned temperature map. We note the follow- 
ing: 



E. Belsole et al.: An XMM-Newton observation of the dynamically active binary cluster A1750 



7 



(alblfcXd) 

Fig. 5. Temperature and entropy maps of A1750. Temperature colour scales are in keV. a: Temperature map obtained 
by extracting spectra in adaptive bins as described in Sect. IT!!^ b: Temperature map derived by applying the wavelet 
algorithm as described in Sect. 14.41 This temperature map was obtained using data from the MOS cameras only, c: 
Specific entropy map. Scale is in arbitrary units, d: Colour coded map of the error in temperature. Pixels with errors 
greater that 1.5 keV have not been taken into account. 



— The hot rectangle to the South of A1750 C in the adap- 
tively binned map is not seen in the wavelet map, but 
this region is at the periphery of the cluster emission 
and has large errors (kr=5.7l}'7 keV). 

— The hot regions appearing in the wavelet map to the 
south and south-west of the centre of A1750 C appear 
to be due to point sources. 

— The hot region between the clusters is recovered in the 
wavelet map, but is compromised by a point source. 
Because of this, we are unable to determine if the struc- 
ture is continuous or clumpy. 

— The temperature structure the the east of the centre 
of A1750 C (discussed in more detail later on) is ro- 
bust, as there are no associated point sources. It is re- 
solved in the wavelet map into two arc-like structures 
of higher temperature. The first is located at ~ 1.'8 
from the centre of A1750 C to the east. The second 
one starts at around 3. '5 in the same direction. These 
two hot regions are separated by a gap of lower tem- 
perature. 

— Finally, the hot region at the Northwest end of the arch 
between the clusters appears to be stable, appearing 
in both maps, although some of this emission may be 
due to a strong point source. 

The western side of A1750 C appears to be relatively 
isothermal. A1750 N displays a quite uniform temperature 
distribution, with an indication of a higher temperature in 
the centre than that observed with the adaptive binning 
map (but see also Section |HJ). 

Note that the wavelet temperature map has been de- 
rived using only the MOS data, while the adaptive binning 
map uses data from all three cameras. As such, we would 
expect the wavelet map to only have approximately half 
the sensitivity of the adaptively-binned map. This loss of 
sensitivity is crucial in regions of low surface brightness, 
such as the middle region. Despite these caveats, the agree- 
ment between the maps is good. 



5. Entropy 

The multi-scale spectro-imagery algorithm described 
above can be applied to estimate the projected spatial 
distribution of other parameters of the intra-cluster gas, 
provided that the expected parameter and its fluctuation 
are estimated within each resolution element. In this pa- 
per, we applied the algorithm to the EPIC-MOS data to 
compute a crude map of the the gas specific entropy , 



Fig. 6. Adaptively-smoothed MOS image in the energy 
band 0.3-6.0 keV. The lines limit the regions used to ex- 
tract 1 dimensional profiles. The image in in logarithmic 
scale. 



which is convenient to define in terms of observed quanti- 
ties as S" oc T/n^/"^, as a function of the intra-cluster gas 
temperature T, and density n^. 

To do so the averaged cluster emissivity per surface 
unit (N) and temperature (T) are fitted simultaneously 
within the different resolution elements, following the 
same procedure as for computing the temperature map. 
Since the cluster emissivity (N) scales as to the square 
of the density of the intra-cluster gas (n), both quantities 

leads to a crude estimator of the averaged specific entropy, 
Q - ^_ 

>-> — 7V1/3 ■ 

Just like for the temperature map, the 'entropy' map 
is obtained by computing maps of S" ± cr^ at five different 
scales, leading to a wavelet transform. Then the wavelet 
coefRcients are thresholded in order to remove the noise 
contribution, again following a la significance criterion. 
The entropy map, shown in Fig. 



is very similar to 
entropy maps obtained in numerical simulations of clus- 
ter mergers (e.g. Ricker & Sarazin '200T). We observe a 
lower entropy in the centres of the two clusters, with an 
increase towards the external regions. The entropy distri- 
bution of A1750 N is very elongated towards the north, 
mainly following the emissivity distribution. We observe 
that the eastern side of A1750 C displays an enhancement 
in the entropy which does not follow a spherical distribu- 
tion. The projected entropy distribution is quite turbulent 
in this region. The entropy in the middle region is also 
higher than in the centre of the two clusters, as expected 
from shock heating. 



6. One dimensional analysis 

The two dimensional analysis above has revealed a wealth 
of new detail in both the morphology and temperature of 
the gas in this cluster. In this Section we use the two di- 
mensional results to divide each subcluster into discrete 
regions, with the aim of undertaking a classical ID anal- 
ysis for comparison with previous results. The regions we 
have chosen are shown in Figure |H| These regions are: 



^ more precisely, the definition of "entropy" is the logarithm 
of the quantity above plus a constant (see Ponman et al. I2UU3I 
and references therein) 
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— Region f (Rf) corresponds to the part of Af750 N 
which is relaxed-looking and fairly isothermal. Rl has 
been defined by excluding all data between position 
angles 77° — 330° (anticlockwise from N, centred on 
a = 13''3f"00";(5 -0f°46'46" ), related to all hot 
regions connected with the interaction. 

— Region 2 (R2) is that part of A1750 C which is relaxed- 
looking and very isothermal. R2 is defined by a sec- 
tor between position angles 205° — 354°, anticlockwise, 
centred on a = f 3''30™49^88f , S = -01°5f '46'.'70, the 
A1750 C emission peak. 

— Finally, Region 3 (R3) is delimited by a sector between 
position angles 65° — 205° (anticlockwise), centred on 
the A1750 C emission peak. R3 includes the zone ~ 4' 
the east of the centre which shows the hot, arc-like 
temperature feature and isophotal compression. 

6.1. Data analysis 

6.1.1. Surface brightness profiles 

In each Region, we extracted azimuthally averaged sur- 
face brightness profiles from the observation and the cor- 
responding Region in the blank-sky background. Photons 
were binned in circular sector annuli of width 5" centred 
on the peak of the emission. We use the energy band [0.3 - 
1.4] keV to minimise the dependence of the emissivity on 
the temperature; the upper energy limit was set to avoid 
the instrumental Al and Si lines^. 

The background surface brightness profiles were sub- 
tracted from each corresponding cluster profile, using the 
appropriate [10. - 12.] ([12. - 14.]) MOS(pn) keV normali- 
sation factor. Each camera was treated separately. In all 
cases, the MOS and pn profiles agree well barring the ex- 
pected normalisation differences, and so in each case the 
MOS and pn profiles were co-added. 

Except in one case (R3), the surface brightness profile 
was binned such that (i) at least a S/N ratio of 3(t was 
reached, and (ii) the width of the bins increased with ra- 
dius, with A(9) > 0.19. The logarithmic binning ensures 
a roughly constant S/N in the outer parts of the profile. 

6.1.2. Temperature profiles 

Radial temperature profiles were produced for each Region 
using the same mask as was used for the corresponding 
surface brightness profiles. Spectra were extracted in cir- 
cular sector annuli centred on the peak of the emission. 
The widths of the annuli were chosen so that a minimum 
5(T detection was achieved in the [2.0 - 5.0] keV band. 
A minimum width of 30" was also imposed, correspond- 
ing to the 90% encircled energy radius of the MOS PSF. 
The blank-sky background spectra were extracted from 



^ We checked the emissivity-temperature dependence of the 
FeL blend by extracting profiles with the [0.9 - 1.2] keV band 
excluded. The best fitting /3-model parameters did not change 
significantly in this case, thus to maximise the S/N, we use the 
entire [0.3 - 1.4] keV band. 



the same regions and subtracted fom the source spectra. 
The background subtracted spectra were binned to 3a to 
allow the use of Gaussian statistics. 

Spectra from the three cameras were fitted simultane- 
ously with an absorbed MEKAL model, with the column 
density fixed to the best-fit of the corresponding global 
spectrum (Tab.[2Jl. The abundance for the fits was frozen 
at the best-fit global value for whichever cluster was under 
consideration'^. The MOS spectra were fitted in the range 
[0.3-10.0] keV, and the pn spectra in the range [0.4-10.0] 
keV. 

We stress that we have not attempted to undertake a 
full PSF and deprojection analysis and that the results 
presented here are projected temperature profiles. PSF ef- 
fects will be important if the cluster possesses very signif- 
icant temperature gradients, or if the surface brightness 
profile is very peaked. We have taken care to minimise 
PSF effects by ensuring that the bins have a minimum 
width greater than the 90% encircled energy radius of the 
PSF, and we have already seen that the temperature gra- 
dients, while significant, are relatively mild (not nearly so 
strong as in a cooling fiow, for example). 

The temperature plots below are presented with the 
upper X-axis in arc-minutes. The lower X-axis shows 
the radius in terms of the virial radius, calculated us- 
ing the average temperature from the global spectral fit 
(Sect. 14.1(1 . and the r2oo — T relation of Evrard, Metzler 
& Navarro 11996)1 at the cluster redshift, viz: 



r-200 = 3690 (T/IO keYy/^{l + z)-^/^. 



6.1.3. Gas density profiles 



(1) 



The surface brightness profile of each Region was fit- 
ted with various parametric models, all of which were 
convolved with the XMM-Newton PSF (Ghizzardi 2001, 
Griffiths & Saxton 2002) and binned in the same way as 
the observed profile. The surface brightness profile, at low 
energy, is generally well represented by a /3-model: 



s{e) = So (i 



-3/3+0.5 



(2) 



where 5*0 is the central intensity, 9c the core radius and 
P the slope. In the hypothesis of an isothermal cluster, the 
two parameters 9c and (3 are related to the gas density 
profile by: 



nga.,ir) = no ( 1 + ^ 



-3/3/2 



(3) 



In addition to the standard /9-model above, we 
have also used the BB parametric model (Pratt & 
Arnaud F2002|l . which is a double isothermal /3-model which 
assumes that both inner and outer gas density profiles can 

We did fit the radial profiles with spectra where the abun- 
dance was a free parameter: at 90% confidence level however, 
the temperature profiles do not show significant differences 
from those obtained with the abundance frozen at the global 
value. 
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Fig. 7. The surface brightness profile of A1750 N, Region 
1, binned as described in the text. It has been background 
subtracted and corrected for vignetting. The profile is 
shown with the best fit BB model; the residuals are shown 
in the bottom panel. The dash-dotted line shows the best- 
fit single /3-model to the outer regions. 
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Fig. 8. Temperature profile of A1750 N, Region 1. The 
straight line is the best fit global spectrum temperature, 
the dash-dotted lines denote ±10% of this value. Errors 
are at la for one interesting parameter. 

Table 3. Results of the BB model analytical fits to the gas 
surface brightness profiles, errors are 90% confidence. 



Parameter 


Rl 


R2 




A1750 N 


A1750 C W 


/, 1/2 -3n 

'^H,o(h5o cm ) 


7.68 X 10"^ 


5.89 X 10"^ 


rc (' / kpc) 


l!96+^',**V 254.4 

-0. 17 
'-'•>J^-0.04 


2!03+'''/^ / 263.5 

-0.12 

69+° °^ 


13 


Rent (' / kpc) 

rc,in (' / kpc) 


2/Q4+2!64 , g 

-0.23 

0:i2_o.o4" / 15.6 


I'OSlaie / 136.3 
0:i0_o.o2" / 13.0 


xVd.o.f 


57.4/35 


42.2/42 



The maximum value of Vc in is fixed to 1'. 



each be described by a /3-model, but with different pa- 
rameters. The boundary between the two regions is a free 
parameter of the fit, and the density distribution, and its 
gradient, are continuous across the boundary. 

6.2. A 1750 N; Rl 

A single /3-model as described above, while being a rela- 
tively good description of the outer regions, is not a good 
description of the entire profile (x^ = 122.5 for 37 degrees 
of freedom (d.o.f)). As is often found in these cases, pro- 
gressively excluding the central regions improves the fit. 
The best fit /3-model, obtained by excluding the inner ~ 1' 
yields — 48.9 for 27 d.o.f.; in this case /3 = 0.46 and 
Tc = 1.'4. We thus used the double isothermal BB model 
to fit the profile. This resulted in a much better fit to the 
entire radial range: x^ = 57.4 for 35 d.o.f. The BB model 
fit is shown in Fig.d best- fit results are shown in Table|31 
The best-fit outer /3 value of the BB model fit is in 
excellent agreement with the 2D analysis in Sect. 13.21 (re- 



sults in Table and with the ROSAT-derived results 
presented by Donnelly et al. H2001|l . The fact that this is 
so, when we are not fitting exactly the same regions, is an 
indication that the results derived here are very robust. 

The radial temperature profile is shown in Fig.|S| note 
that this is the first such profile for this cluster. The clus- 
ter is detected with excellent S/N up to ^ 0.4r2oo, the 
limit being due to the edge of the XMM-Newton FoV. 
The temperature profile in this Region shows a surpris- 
ing amount of variation. However, while the variations are 
significant in terms of their errors, they are actually quite 
small. For instance, the lowest temperature annular sec- 
tor has kr = 2.6 keV, and the highest kT = 3.3 keV, an 
absolute variation of 0.7 keV. 

There is a dip towards the centre, which may be due 
to a cooling flow (discussed further below in Sect. 17. 3|) . 
From the second to the final annular sector, there is a 
smooth, gentle decline broken only by the fifth annular 
sector. However, the general features of this temperature 
profile are in good agreement with the spatially-resolved 
temperature map (Sect. ^3 . 
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Fig. 9. The surface brightness (top) and temperature (bottom) profiles of Regions 2 and 3 of A1750 C. The right hand 
panel in each case shows R2, the left hand panel shows R3. For R2, the dot-dash line if the best-fit /3- model and the 
solid line is the best-fit BB model. The residuals are also shown. The dotted line overplotted on R3 is the best-fit BB 
model for R2. 
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6.3. A 1750 C; R2 

R2 extends to the eastern side of A1750 C, and cor- 
responds to the part of the cluster which shows the 
least amount of structure in isophotes and appears quite 
isothermal in the temperature map. 

Once again, a single /3-model fails fully to describe the 
data (x^ — 56.7 for 44 d.o.f.), although the central excess 
is far less evident than for Rl. The best fitting single f]- 
model requires the exclusion only of the central 0.'7, and 
yields = 36.6 for 37 d.o.f. A BB model with parameters 
given in Table |3| is an excellent fit to the entire profile 
(X^ = 42.2/42 d.o.f). The profile, together with the best 
fitting BB model, is shown in the right hand top panel of 
Fig. El 

The best-fit outer /3 value of the BB model fit is in 
good agreement with the 2D analysis in Sect 13.21 and 
with the ROSAT-derived results presented by Donnelly 
et al. (|2()()l|l . Once again, this gives us good confidence in 
our results. 

The temperature profile for R2 is shown in the right 
hand bottom panel of Fig. EJ this is the first temperature 
profile for this cluster. 

The cluster is detected with excellent S/N up to ^ 
0.6r2oo- In this Region, as expected from the temperature 
map, the temperature profile is quite smooth, and shows 
little variability. 

6.4. A1750 C; R3 

Here we come to the region of A1750 C which displays 
isophotal compression and, compared to R2, considerable 
temperature structure. 

The surface brightness profile of R3 is shown compared 
to the profile from R2 in the left-hand top panel of Fig.O 
in order to avoid smoothing of the discontinuity, we have 
not adopted logarithmic binning.. The difference between 
the profiles is striking: the R3 profile is more peaked, as is 
shown by the overplotted best-fitting BB model to R2. In 
addition, there appears to be a change of slope between 
3 and 4 arc-minutes from the centre, the shape of which 
may indicate a discontinuity in the gas density profile. 

For simplicity, we assume spherical symmetry. To 
quantify the discontinuity, we fitted the surface bright- 
ness profile of R3 with a radial density model composed 
of a /3-model and a power law separated by a jump: 



n(r) 



where 



n2ir) = ^2(1 + 



-3/3/2 



ni(r) = Al H- 



if r < Tc 



if r > 



(4) 



nil 



= AT: 



jump 



In this model, rt2 and ni are the gas densities before 
and after the discontinuity, respectively. The free param- 
eters of the fit were the slope /? and coreradius of the f3- 
model, the slope a of the power law, and the position (rcut) 



A 1750 




3. 4,5 
radius (arcmin) 



0.1 



/S model + powerlaw 



Radius (arcmin) 

Fig. 10. Top: Surface brightness profile of A1750 C, 
R3: data points and density model used to fit the data. 
Bottom: Model used to fit the surface brightness profile 
(see Eqn. 2J|. The model is normalised to 1 at the centre. 



and amplitude Xjump of the density jump. The profile was 
fitted between bin 2 (10") and bin 60 (6'). The best fit den- 
sity model is shown both alone and superimposed on the 
data in Fig.^| The goodness of the fit gives an excellent 
51.4/53 d.o.f. The best fit parameters are: (3= 0.41, 
rc=0.'6, q;=0.85. The jump position is rcut = 3.44' and the 
amplitude of the density jump is A'jump = 0.78^^1^, i-e., 
a 20% jump in density. 

We also changed the bin range within which we fit the 
model and also used a second model composed of 2 power 
laws. In both cases we obtain consistent results within 
the errors, including jump amplitudes, indicating that the 
results are robust. 

The orientation of the discontinuity with respect to 
the line of sight is unknown: projection effects may act to 
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reduce the apparent density jump at the discontinuity. If 
we do not see the discontinuity exactly perpendicular to 
the line of sight, the sharpness of the edge will be reduced. 
In other words, it is likely that our derived value for the 
density jump is a lower limit. 

The radial temperature profile for R3 is shown in the 
bottom left panel of Fig.^ Here, the ring edges were cho- 
sen in order to match the density jump described above, 
and to have at least a 5cr signal above the background in 
the [2.0 - 5.0] keV band. 

The temperature profile for this region is quite vari- 
able. Even bearing in mind that both temperature pro- 
files in Fig. 1^1 are plotted with a logarithmic x-axis, it is 
obvious that the profile of R3 is more variable than that 
of R2. The maximum of the temperature profile at lib 
(in fact a la deviation) corresponds to the peak of the 
inner high-temperature arc seen in the temperature map 



Radius (orcmln) 



(Figure 5(b) I. This temperature peak does not correspond 
to any significant structure in the surface brightness pro- 
file, even if the surface brightness profile at this distance 
from the centre does display some level of discontinuity. 
It is worth noting that within the errors we do not detect 
any temperature discontinuity at 3'. 4, the position of the 
jump in the density profile. We further note that there is 
a significant temperature jump between annuli 5 and 6. 
However this does not correspond to any significant jump 
in the observed surface brightness profile. 

6.5. Entropy profiles 

With a view to making a comparison with relaxed clus- 
ters, we determined the entropy profiles of Regions Rl, R2 
and R3. The entropy was determined from the BB analyt- 
ical model fit to each gas density profile and the observed 
(projected) temperature profiles, taking, as is now cus- 
tomary in this field, S = Tjr^J'^ . For R3, we use the gas 
density model described above in Sect. 16.41 limiting the ra- 
dial values to the fit limit of r < 6'. The resulting entropy 
profiles are shown, plotted in terms of the virial radius, in 
Fig. ^2 Typical errors, corresponding to the error in each 
bin of the temperature profile, are indicated. The over- 
plotted line shows the S oc T^^^ behaviour expected from 
shock heating (Tozzi & Norman I2001|l 

7. Discussion 

7.1. The current merging event 

In this Section we discuss a possible scenario for the cur- 
rent merging event (i.e., that between A1750 N and A1750 
C) based on simple assumptions. 

The weak enhancement in temperature between the 
two clusters (Fig. 5(a) I is a strong indication that the 



clusters are just beginning to interact. 

We can attempt to calculate the Mach number and 
the collision velocity from the observed temperature vari- 
ations. In this analysis, we follow the simple assumptions 
outlined in Markevitch et al. ((1999.) . viz: that the merger 
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Fig. 11. The entropy profiles of Regions Rf , R2 and R3. 
The dashed line corresponds to S" ct r^-^, with an arbitrary 
normalisation. The profile for R3 has been limited to the 
range of our fit to the surface brightness (r < 6'). 



is symmetric in the plane of the sky, that the shocked 
gas is on average at rest with respect to the centre of 
mass, and that the shocked gas is nearly isothermal and 
in equipartition. The Mach number and collision velocity 
can then be obtained by applying the Rankine-Hugoniot 
jump conditions for a one-dimensional shock, following 
Landau & Lifshitz 1)19591 see also Markevitch et al. 119991 

Sarazin mm-. 



1 

c ~ 



Ti 



3C 



T2 

Ti 



1/2 



Ti 



(5) 



(6) 



where Eqn. [3 corresponds to Equation 2 of Markevitch 
et al. 1)1999(1 when 7 = 5/3; Ti, T2 are the pre-shock and 
post-shock temperatures and C — n2/ni is the shock com- 
pression. 

For this rough estimate, we assume that the pre-shock 
temperature (kT = 3.10 ±0.14 keV) can be approximated 
by the average temperature on the sides of each cluster 
symmetrically opposite the shock region (see Markevitch 
et al. 11999(1 . and the post-shock temperature (kT = 5.12± 
0.73 keV) is given by the observed global temperature of 
the middle region (Sect. ll!T|l . We obtain C = 1.89 and a 
Mach number of = 1.64. 

Assuming that the surface brightness scales as the 
square of the density, using the above results we expect 
an increase of ~ 3.6 in the surface brightness between the 
clusters. This is not a huge effect. To compare directly 
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with the observations we integrated the surface bright- 
ness perpendicular to the axis joining the cluster centres, 
and compared this to the same integration of an image of 
the best-fitting BB models. We then computed the ratio 
in brightness between the middle region and the region 
symmetrically opposite the clusters (in a similar fashion 
to the determination of the Mach number above) , for both 
the image and the model^. The range of ratios we found 
(1.5 - 6.2) brackets the expected increase. 

To first order, the time until first core passage (tmorgc) 
between A1750 C and A1750 N can be obtained by consid- 
ering the distance of the two components and the distribu- 
tion of the temperature and the density, even if projection 
effects are the major limitation in the accuracy of such an 
analysis. 

The centres of A1750 C and A1750 N he at a pro- 
jected distance of 675 h'^Q kpc. If we assume that the 
difference in radial velocity is exclusively due to the infall 
of A1750 N onto A1750 C (in the reference frame of A1750 
C), the two systems would be so close — well within the 
virial radius of each cluster — that we should observe far 
larger density/surface brightness distortions and temper- 
ature enhancements in the middle region (as seen in e.g., 
Fig. 3 of Ricker & Sarazin 2001). On the other hand, if 
we assume that the physical distance is determined by the 
Hubble flow, meaning that the motion of the two clusters 
is exclusively radial, the difference in their mean radial 
velocities would yield a distance between the two clus- 
ters of ~ 20 h^Q Mpc. This corresponds to nearly 10 
times the virial radius of cither cluster and thus the two 
should not display any sign of interaction whatsoever, in 
disagreement with the results described in detail above. 
The real distance surely lies in between these two values, 
and judging by the observed characteristics, is likely to be 
less than the virial radius of A1750 N (~ 1.7 Mpc), using 
the r2oo - T relation of Evrard et al. I|1H96I - Eqn.IU). 

Using Equation 1 of Markevitch et al. (1999), we cal- 
culate a collision velocity of 1400 km s"^, which is roughly 
of the same order as the radial velocity difference between 
A1750 C and A1750 N. For a rough, order of magnitude 
estimate of imcrgc, we first assume the distance between 
the clusters to be the projected distance. In this case, we 
find Emerge = 0.47 Gyr. If we now assume that the distance 
between the clusters is given by the estimated virial radius 
of A1750 N, we find Uergc = 1-3 Gyr. 

We can also compare our results with simulations of 
merger events of systems of nearly equal mass (Roettiger 
et al. 1996, [TM7I Ricker & Sarazin _2Q01i Ritchie & 
Thomas |2()()2I Teyssier 2002), which resemble very closely 
the observed characteristics. We thus expect the two clus- 
ters to lie at a real distance close to, but somewhat lower 
than, their virial radii, around 1 h^Q Mpc. This would im- 
ply that the clusters will reach core passage some time in 
the next - 0.7 Gyr. 



7.2. A1750 C: an unrelaxed cluster 

Four pieces of evidence lead us to conclude that A1750 C 
is not a relaxed cluster, and further, that the perturbed 
state is intrinsic to the cluster itself and is not connected 
to the nascent merger with A1750 N. These are: 

1. The observed discontinuity in the gas density 
profile. 

A shock front is the most natural interpretation for 
the density discontinuity observed in A1750 C, R3 
(Sect. I6.4|l . In addition, the entropy map and profile 
of R3 (Fig. Ill|l show exactly the type of behaviour 
expected from a shock (the denser side has higher en- 
tropy). 

Heating shocks are expected on theoretical grounds 
when clusters collide. Furthermore, such shocks are 
expected to be relatively weak, with a distribution 
of Mach numbers peaked at X = 1.4 (Gabici EESI)- 
Shocks were proposed to explain the peculiar tem- 
perature and density structures observed in Cygnus 
A and A3667 (Markevitch et al. 1999) with ASCA 
and ROSAT. A clear detection of a shock wave was 
obtained with Chandra in the IE 0657-56 cluster 
(Markevitch et al. "2002). Other weaker shocks were 
observed in A3667 (Vikhlinin, Markevitch & Murray 
nUUT^ and A85 (Kempner, Sarazin & Ricker I^IMjl . 
By applying the Rankine-Hugoniot jump conditions 
(Eqn.EJ we can obtain the Mach number of the shock 
in A1750 C with respect to the ambient gas (assum- 
ing again an adiabatic index for monoatomic gas 5/3). 
These conditions are valid if the gas within the shocked 
region is nearly isothermal and equipartition between 
electrons and ions applies. From the density jump in 
Sect.EHwe obtain M = l.ldtom- The expected tem- 
perature ratio for a shock with such a Mach number is 
T2/T1 = 1.18. Comparing the pre-shock ri=4.28l°;^^ 
keV and post-shock T2 = 4.5915]^^ keV temperatures 
(see Sect.ES, we obtain a ratio Tz/Ti = 1.07(±0.21)5, 
and thus, within the errors, the two independent mea- 
sures are consistent. 

From the jump conditions we can also infer the pres- 
sure before and after the shock, obtaining a pressure 
jump of a factor 1.5. This demonstrates that the gas 
in cluster C is not in pressure balance and thus hydro- 
static equilibrium conditions should not be applied. 
The geometry of the shock suggests that it is related 
to a merging event that A1750 C has suffered in the 
past. The merger velocity is Vmerge ~ 500 km s~^, 
and thus taking the distance from the centre, a rough 
estimate of the time of the A1750 C merger gives 
tmerge ~ 0.9 Gyr. By comparing with simulations, the 
thermo-dynamical status this corresponds to an ad- 
vanced phase of merger, at least 2 Gyr after core pas- 
sage (Ricker & Sarazin 1200 Ifl . which also corresponds 



* The definition of the outer regions is somewhat arbitrary 
since we do not know the real distance between the clusters 



^ Assuming symmetrical errors for Tl and T2 



14 



E. Belsole et al.: An XMM-Newton observation of the dynamically active binary cluster A1750 



to a sound crossing time^. Given our approximations, 
possible projection effects and the physical changes in 
the medium during the merger, and the difficulty of 
comparing with simulations, we cannot really be more 
precise. 

The observed temperature jump between bins 5 and 
6 in the R3 temperature profile brings to mind a cold 
front. However, the propagation of a cold core into a 
hotter medium should also produce a density jump, as 
observed in e.g., A2142 (Markevitch et al. I^MIjl . We 
have not been able to find such a discontinuity at this 
position, even by using profiles extracted in elliptical 
rings. However, given possible projection effects, we 
cannot entirely rule out the possibility of a cold front. 
Obviously, the better spatial resolution of a Chandra 
observation would enable us to provide a definitive an- 
swer to this question. 

The elliptical isophotes at small scale, the offset 
of the core from the expected position given the 
larger scale structure, and the shift of the X- 
ray centroid from the position of the brightest 
cluster galaxy (BCG) . 

These are particularly evident after subtraction of a 
2D /5-model (Fig.HJ, and in Fig. 



Radius (orcmln) 



1(b) 



In their simulations of merging clusters Roettiger et 
al. (1997) predict a late phase of the interaction, af- 
ter core passage, when the dark matter and the gas 
separate and the gas sloshes in a potential which is 
starting to re-establish a new equilibrium configura- 
tion. This phase is long in the timeframe of the merger 
event (at least up to 5 Gyr in the Roettiger et al. simu- 
lations) . This sloshing of the gas can produce an elon- 
gated distribution as observed in A1750 C as well as a 
centroid shift (sec also Roettiger et al. 1996). A simi- 
lar gas distribution is observed in A1795 (Markevitch 
et al. 2001; Fabian et al. 2001; Ettori et al. 2002). In 
this case, however, the extent of the sloshing gas is 
1/2 of that of A1750 C and is more closely connected 
with the central galaxy and the cooling flow. Note that 
Markevitch et al. (2001) and Ettori et al. (2002) con- 
cluded that the centre of A1795 is not in hydrostatic 
equilibrium even at larger scale. This is thus further 
strong observational evidence that A1750 C is in the 
late phase of an old merger event. 
The lack of evidence for a cooling flow. 
The flat central temperature gradient (Fig. ^ argues 
quite strongly against the presence of a cooling flow. 
We confirm this by computing the cooling time us- 
ing Eqn. □ This gives tcooi = 1-05 x 10^° yr, of the 
same order as the Hubble time at the redshift of the 
cluster. Simulations of merging clusters suggest that 
only nearly head-on major mergers between objects of 
nearly the same mass can destroy cooling flows (e.g., 
Ritchie & Thomas 2002). Thus the lack of cool gas in 



^ ts= 6.610** {T/10^K)-^/^{D/Mpc) 
diameter (see Sarazin Fiybtil ea. 5.54) 
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Fig. 12. The scaled entropy profiles of Regions R2 and 
R3 compared with the entropy profiles of A1413 and 
A1983 (Pratt & Arnaud 2002, 2003). The profiles have 
been scaled with using the empirically-derived relation 



S oc T^' 



1-0.65 q£ pQnjj-^aji^ Sanderson & Finoguenov H2003|l . 
The dashed line corresponds to oc r^'^, with an arbitrary 
normalisation. 



where D is the cluster 



the centre of A1750 C is likely the result of an old - 
and rather violent - merger event. 
The entropy profile, which appears to diffter sig- 
nificantly from other high quality XMM-Newton 
profiles of relaxed clusters. 

We show in Fig. El the entropy profiles of R2 and 
R3, compared with the profiles of the relaxed clus- 
ters A1983 (kT = 2.1 keV; Pratt & Arnaud 2003) and 
A1413 (kT = 6.9 keV; Pratt & Arnaud 2002). The 
profiles have been scaled with the empirically-derived 
scaling relation {S oc T^°-^^) of Ponman et al. PMji . 
While both R2 and R3 appear approximately to con- 
verge to the S oc r^'^ value expected from shock heat- 
ing (Tozzi & Norman I2001|l at a radius of ~ 0.2r2oo, 
each exhibits a very high scaled central entropy value 
for its temperature (or mass), compared to the scaled 
entropy of the relaxed clusters. 

In form, the entropy profiles of R2 and R3 are qual- 
itatively similar to those found in the merger simula- 
tions of Ricker & Sarazin H2001|l . The entropy map is 
also reminiscent of their simulations. The increase in 
core entropy due to a merger has also been seen in the 
simulations of Ritchie & Thomas l|2002|l . The increase 
appears to be the result of the mixing of shocked, high 
entropy gas from the outer regions into the low en- 
tropy core of the cluster. Thus the entropy profile is 
another strong piece of evidence that A1750 C has not 
yet relaxed from a previous merger event. 
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Crucially, all of the evidence listed above solely con- 
cerns A1750 C. These data strongly suggest that the clus- 
ter is in the late phase of a merger for which the core pas- 
sage happened 1-2 Gyr ago. The cluster is in the process 
of re-establishing equilibrium, but this is likely to be in- 
terrupted by the new merging event occurring with A1750 
N. 

7.3. A1750 N 

The morphological analysis of A1750 N (Sect. I3.2|l has 
shown the elongated shape of the cluster, as well as an ex- 
cess above a standard /3-model which appears considerably 
off-centre with respect to the overall cluster morphology. 

In discussing the form of this cluster we should first 
mention that at this off-axis angle PSF effects are expected 
to become important. In order to have a qualitative idea 
of this effect, we have simulated the whole observation 
in SCISIM, using as input the best-fitting ID /3-model 
values listed in Table|21 The resulting image confirms that 
emission from A1750 N should indeed be smeared by the 
PSF, preferentially in the azimuthal direction. 

With this in mind, we now turn to the ROSAT image 
of the cluster, shown in Fig. 1131 and discussed in more 
detail below (Sect.[73Il- It is evident that A1750 N seems 
elongated in the NE direction; however, given the larger 
field of view of ROSAT, we notice that there is another 
extended source just to the North of the cluster. We can 
thus speculate that there is some contamination by emis- 
sion from this source, and indeed the 2D /3-model analysis 
of the XMM-Newton observation suggests extended resid- 
uals in the direction of this source. Again we caution that 
PSF effects play a role at this distance from the centre. 

Finally, we cannot entirely rule out that A1750 N has 
been disturbed, or is being disturbed, by interactions. It 
is possible that the observed isophotal compression to the 
South of the centre (see figure is due to the interaction 
with A1750 C (suggesting a non-zero impact parameter), 
but with the present observation we cannot definitively 
say whether there is interaction with the source to the 
North. 

Donnelly et al. (2001) speculate on the presence of a 
cooling flow in A1750 N. This is qualitatively in agreement 
with what we observe: a temperature drop is observed in 
the core of A1750 N (Fig.lSJ, together with an excess above 
a standard /3-model. 

In order to verify the presence of a possible cooling flow 
in A1750 N, we calculated the cooling time using (Sarazin 
1986): 

-°'V(^)"(^)-. (7, 

Using the central gas density obtained in Sect. IS. 21 
ne(0)=7.68 X 10-3 cm'^; Eqn.|7| gives tcooi = 6.4 x 10*^ 
yr, about half of the age of the Universe at the redshift 
of the cluster. This is consistent with a weak cooling flow. 
We thus fitted the MOS spectrum of the central bin with 
more complicated models: 



Table 4. Multi-temperature and CF fits of the inner annu- 
lus of A1750 N. The _F-test is computed against the fit for a 
single temperature absorbed MEKAL model. Errors are 90% 
confidence. 



Parameter IT 2T CF 



kTi (keV) 


9 77+U.45 
^- ' ' -0.35 


o Q1 +U.V9 


o 90 + 1.4 
'J-^O_0.46 


kTz (keV) 




' •J-0.23 




Z/Zq 


n 9+0.22 
'-'•^-0.16 


0.251°:?^ 




M {Modot yr 
xVd.o.f 






4.97 (< 9.6) 


108.22/100 


105.15/98 


105.77/98 


Fproh 




75.6% 


67.4% 



The low temperature of the CF model Tiom was 0.71 keV. 



— The sum of two MEKAL models absorbed by a col- 
umn density fixed to the global spectrum fit value. The 
abundances of the two components are tied together. 

— The sum of a MEKAL model and a cooling flow model 
(MKCFLOW), with flxed absorption and abundances 
of the CF model tied to the thermal model. Also the 
upper temperature of the CF model is tied to the tem- 
perature of the MEKAL model. 

The results are shown in Table 0] together with the 
simple MEKAL model results for this bin. There is 
not a significant probability that the 2T model or the 
MEKAL-I-CFLOW model are a better representation of 
the data. All fits are very good and the small mass depo- 
sition rate found by fitting the CF model is in agreement 
with a small cooling flow. 

Apart from the cooler core region, A1750 N displays a 
fairly uniform temperature distribution. With the present 
data, all evidence suggests that the only physical process 
which is disturbing A1750 N today is the merger with 
A1750 C. However, we stress that a conclusive analysis 
needs an observation of the whole cluster, and the data 
presented here are limited by the off-axis position in the 
XMM-Newton field of view (as discussed further below in 
Sect. mji . 

7.4. At larger scale 

The XMM-Newton FoV has allowed us to observe the two 
main clusters of A1750 in one pointing. However, other 
structures have been observed at larger scale with Einstein 
(Jones & Forman I1HH9(I . Figure El shows a smoothed 
ROSAT/PSPC image obtained from the ROSAT archive: 
the superimposed circle indicates the XMM-Newton FoV. 
To the south of A1750 C we observe another obviously 
extended source. This source has been observed in the op- 
tical and is discussed in Beers et al. poyijl . where it is 
called A1750 S; it lies at a projected distance of ~ 2 Mpc 
from A1750 C and it is at the same redshift as A1750 C. 



Fig. 13. ROSAT/PSPC image smoothed with a Gaussian 
filter of (T = 20". The white circle traces the limit of the 
XMM field of view. 
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Beers et al. H1991|l . using optical virial mass estimators, 
suggest that A1750 S contributes ~ 2% to the total mass 
of the system. It is worth noting that the gas distribution 
to the south-western side of A1750 C appears elongated 
in the direction of A1750 S. 

Fig. ^1 also clearly indicates a further extended source 
just to the north of A1750 N, which is outside the FoV of 
our XMM-Newton observation. We recall that the subtrac- 
tion of a 2D /3- model (Sect. suggested that there are 
further residual structures to the north of A1750 N, likely 
linked to this source. Other possibly extended sources are 
visible further to the north of A1750 N in the PSPC im- 
age. All these sources are aligned in the direction of the 
line joining the centres of the 3 main clusters, strongly 
suggesting a filamentary structure. Einasto et al. H2(JU1|I 
found that A1750 belongs to a rich super cluster com- 
posed of 7 members. We searched clusters in the redshift 
range 0.08-0.09 in a 16° x 16° field centred on the centre of 
this super cluster and found 13 clusters, including A1750 
which is located at the north-eastern edge. 

Given this observed large scale structure, it is not sur- 
prising that A1750 C has suffered one (or several) previous 
mergers along the putative filament and indeed, statistical 
studies suggest a correlation between the level of substruc- 
tures and the environment density (Schuecker et al. 12001(1 . 
This adds a further piece of evidence in favour of our in- 
terpretation. This may well be the case for A1750 N as 
well, where its elongated gas distribution may also be the 
result of accretion along the filament. 

7.5. Caveat emptor 

The scenario outlined above explains in a self-consistent 
manner the data we have analysed. In particular, the 
observed temperature and surface brightness structure 
would be difficult to explain if the clusters have already 
past one another (as suggested for A3528 by Gastaldello 
et al. 12003(1 . We stress, however, that projection effects can 
be very important (e.g., as shown by Roettiger et al. 1997), 
and comparison with numerical simulations is necessarily 
a qualitative exercise. For a better understanding of this 
complex system, two more XMM-Newton pointings, one 
each centred on A1750 N and A1750 S, would be very 
useful, as would deep optical observations. 

8. Summary and Conclusions 

We have reported the XMM-Newton GT observation of 
the merging cluster A1750. The main conclusions of this 
work may be summarised as follows: 

— In the morphological analysis with 2D /3-models, we 
detect excess and off-centre emission in the cores of 
both A1750 N and C. We do not find any significant 
substructure in the region between the two clusters 
after subtraction of the /3-models. However, on close 
examination of the surface brightness in the region be- 
tween the two clusters, we detect an increase of the 



same order as expected if the region between the clus- 
ters is a weak shock region (Sect. 17.1(1 . We further de- 
tect excess residual emission to the north of A1750 
N, which is likely (from comparison the XMM-Newton 
and PSPC images) related to the larger scale structure 
of the system. 

— We have produced temperature maps by applying two 
different algorithms. We measured a temperature in- 
crease of order 30% in the region between A1750 N and 
C. A1750 N exhibits a relatively smooth, uniform tem- 
perature distribution, but there are significant temper- 
ature variations within A1750 C, which appear intrin- 
sic to A1750 C and not connected to the merger be- 
tween the two clusters. 

— We use the 2D information to select interesting regions 
for a ID analysis. For the regions corresponding to the 
relaxed, isothermal-looking parts of each cluster (Rl 
and R2), the ID and 2D /3- model parameters are in 
excellent agreement with each other and with previous 
ROSAT analysis, although the XMM-Newton data re- 
quire an additional central component properly to fit 
the profiles. Both temperature profiles appear to de- 
cline modestly (< 10%) with radius. 

— The density profile of the disturbed side of A1750 C 
(R3) is better described by a /?-model and a power 
law, separated by a jump of order 20%. The tempera- 
ture profile of this region is very irregular, but appears 
constant across the jump. 

— Compared to the scaled entropy profiles of relaxed 
clusters, the scaled entropy profiles of both regions of 
A1750 C exhibit a high scaled entropy value in the 
centre, similar to what is observed in merging cluster 
simulations. Only beyond ^ 0.2r2oo docs the profile of 
R2 attain the S oc r^'^ expected from shock heating. 

— In the assumption of an head-on merger between 
A1750 C and N we calculated the Mach number and 
velocity of the merger from the measured temper- 
atures. With the resulting values (A^ = 1.64 and 
Vmerg ^^-"1400 km s~^), wc estimate that the two clus- 
ters will reach core passage in — roughly — less than 
1 Gyr. 

— By applying the Rankine-Hugoniot jump condition to 
the measured density jump in A1750 C (R3) we con- 
firm the detection of a weak shock of M = 1.2, in 
agreement with the measured temperature of the pre- 
shock and post-shock gas. Under simplified assump- 
tions and by comparison with numerical simulations, 
we estimate that a merger — intrinsic to A1750 C — 
has occurred sometime in the past 1-2 Gyrs. 

— An inspection of the larger scale X-ray image suggests 
that A1750 N and C are the two main clusters within a 
filament which includes at least another small cluster 
at the same redshift of A1750 C (A1750 S). At still 
larger scale, A1750 lies in a filamentary supercluster 
containing, in total, 13 clusters. 

We thus confirm that A1750 is a merging cluster in an 
early phase, when the two units have just started to inter- 
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act and stand at a distance slightly lower than their virial 
radii. However, the global dynamical status of the system 
is far more complicated than expected. We conclude that 
the main cluster (A1750 C) has already suffered a pre- 
vious merger and is now in the phase of re-establishing 
equilibrium. This phase will be interrupted by the cur- 
rent merger with A1750 N. As far as we know, this is the 
first time such complex dynamical signatures have been 
observed in clusters at this merger phase. 

The present day morphology of clusters may thus de- 
pend not only on on-going mergers (or the last merger) 
but also on the more ancient merging history, especially 
in dense environments. This has to be taken into account 
in the interpretation of statistical studies of cluster mor- 
phology. 
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